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Abstract— The paper presents a certainty equivalence state 
feedback indirect adaptive control design method for the sys- 
tems of any relative degree with unmatched uncertainties. The 
approach is based on the parameter identification (estimation) 
model, which is completely separated from the control design 
and is capable of producing parameter estimates as fast as 
the computing power allows without generating high frequency 
oscillations. It is shown that the system’s input and output 
tracking errors can be systematically decreased by the proper 
choice of the design parameters. 


I. INTRODUCTION 


Control design for systems with unmatched uncertainties 
is a challenging task, and the main design method is based 
on the backstepping technique outlined in [4]. There, it has 
been shown that certainty equivalence principle leads to over 
parametrization, which can be avoided by the departure of 
the certainty equivalence. In this case, the adaptive laws 
enter into the control law, which can result in the high 
magnitude control signals in the case of large adaptive 
rates (fast adaptation). An alternative certainty equivalence 
control design method is presented in [1], which avoids 
over parametrization for linear systems, but not for nonlinear 
systems with the relative degree greater than two. 

In this paper, we present a certainty equivalence indirect 
adaptive control approach without over parametrization for 
the nonlinear systems of any relative degree, which is the 
main contribution of the paper. The approach is based on 
the identification scheme, which is completely separated 
from the control design. To enable a fast adaptation without 
generating high frequency oscillations in the adaptive signals, 
it employs an error feedback term, like in the modified 
reference model MRAC (M-MRAC) architecture introduced 
in [5]. It is shown that the state prediction error converges 
to zero independent of the control design. Moreover, it 
is shown that transient of the state prediction error and 
the combined parameter estimation error can be regulated 
by the proper choice of the error feedback gain and the 
adaptation rate. The control design follows the command 
filtered backstepping procedure [2]. It is shown that the 
input tracking error (difference between ideal control and 
command filtered certainty equivalence control signal) and 
output tracking error can be regulated by the proper choice 
of design parameters. 


Vahram Stepanyan is with Mission Critical Technologies Inc. at 


NASA Ames Research Center, Moffett Field, CA 94035, email: 
vahram.stepanyan @ nasa.gov 

Kalmanje Krishnakumar is with Intelligent Systems Division, 
NASA Ames Research Center, Moffett Field, CA 94035, email: 


kalmanje.krishnakumar @ nasa.gov 


The rest of the paper is organized as follows. In Section 
II, we give the problem statement and the assumptions. In 
Section II, we introduce the identification model and give 
its properties. The control design is presented in Section IV, 
and the controller’s performance is analyzed in Section V. 
A simulation example is presented in Section VI, and some 
concluding remarks are given in Section VII. 


II. PROBLEM STATEMENT 
Consider an uncertain single input single output (SISO) 
system in the parametric strict feedback form [4] (p. 99) 
&1(t) = xo(t) + 0] (21) (1) 
&o(t) = x(t) + 03 Go(a1, 22) 


én(t) = u(t) + 4, %,,() 


with 2(0) = a, where x(t) = [x1,...,¢n]' € R” and 
u(t) € R are the state and the input of the system, 0; € R? 
are vector of unknown constant parameters, and y,; : R’ > 


R”’, «7 = 1,...,n are sufficiently smooth known vector- 
functions. The system (1) is written in the vector form 
z(t) = A(t) + bu(t) + n(t), (2) 


where we denote 


eS On—-1x1 = In-1xn-1 IF fez | On-1%1 


0 O1xn—1 1 


and (t) = [01 1 (21) 3 Po(a1,22)---On Pn(w)]". 

The objective is to design a control signal u(t) such that 
the closed-loop signals are bounded, and the system’s output 
y(t) = x1(t) tracks the output y,(t) = 2,1 of the reference 
model 


&,(t) = A,pwr(t)+b,r(t), x,(0) = x0, (3) 
where A, = A— Bk!, b, = k,b, and r(t) is a piece- 
wise continuous and bounded external command. Here, the 
feedback gain k and feedforward gain k,. are chosen to make 
A, Hurwitz and meet the performance specifications. 

III. IDENTIFICATION MODEL 
In order to generate the necessary parameter estimates, we 
introduce the following identification model 
x(t) = A&(t) + bu(t) + A(t) + cx(t) 
&(0) = £0, (4) 
where &(t) is the state prediction, £(t) = a(t) — &/(t) is 
the state prediction error, c > 0 is a design parameter, 


2 aT ~T ~T 

A(t) = 18; (t)ei(a1) 2 (t)Po(11, 72) --- On HPn(e)]", 
and 6;(t) is the estimate of the unknown parameter 0; for 
each 71 = 1,...,n. These estimates are generated according 
to adaptive laws 


6;(t) = y2;(t)¢;(71,.-.,%%), peeeg , (5) 


where y > 0 is the adaptation rate. The state prediction error 
dynamics do not explicitly depend on the control signal 


a(t) = (A — cl)&(t) + 4(t), (6) 


where we define #(t) = n(t) — A(t) with 6,(t) = 0; — 
6,(t), i =1,...,n being the parameter estimation error. 
Lemma 3.1: The error signals (t) and 6,(t), i = 
1,...,n are globally uniformly bounded, and z(t) > 0 as 
t> ow. 
Proof: Consider a candidate Lyapunov function 


i=1 


1S Toe aes 
V== i (t) + —0, (t)6;(t 7 
O=FXRO+ TOO], -@ 
the derivative of which is computed along the trajectories of 
the prediction error dynamics (6) and the adaptive laws (5) 


V(t) =@" ()(A— D)@(t) 
SPB) lateve, Site) ao 


Substituting the adaptive laws and completing the squares 
results in 


i=1 i=1 

With c > 1, the LaSalle-Yoshizawa theorem ( [4], p.24) 

guarantees that &(t), @;(t), i = 1,...,n are globally 

uniformly bounded, and %(t) — 0 as t > oo. In particular, 

there exists 3; > 0 such that 57", ||6;(t)||? < 6?. a 
Lemma 3.2: If &) = Zo, then the state prediction error 

«(t) satisfies the bound 


@@llc. < Bry”. (9) 
Proof: The proof follows from the fact that ||a&(t)||? < 
2V(t) < 2V(0) < G7 for all t > 0. a 


The next lemma gives the bound on the state prediction 
error when #% # Zo. 
Lemma 3.3: If #9 # Xo, then X(t) satisfies the bound 


#(2)|| < Bye 4 BE 


; (10) 
vy 
where $2 = ,/|2V(0) — Bi), 
Proof: The inequality (8) can be written as 
: —1 
V(t) < -2(c— 1)V(¢) + —=8?, any) 
Y 


which implies that 


ae 


i —2(c—1)t BF —2kt , © 


Recalling that ||z(t)||? < 2V(t), we obtain 


oe | “| —2(c-1)t Bi 
[|@(e)I| < yf |2V(0) — —] en%e-ME+ —, (13) 
of ay 
The bound (10) follows from the fact that a+b < JVa+vb 
for any a>0, b> 0. | 

Since the effect of the initialization of the state prediction 
error decays exponentially with the rate c— 1, and c will be 
assigned to large values in order to damp the high frequency 
oscillations in adaptive estimates for large adaptation rates, 
we set Zo = Xo in the following derivations. 

Next lemma gives the bounds on the error signal 7(t) and 
a tighter bound on «(t), when a(t) and u(t) are bounded 
(which will be provided with the control design in the next 
section). 

Lemma 3.4: Let the estimates @(¢) and 7)(t) be generated 
by the identification model (4) and (5). In addition, let 
a(t) and u(t) be bounded. Then 7(¢) and &(t) satisfy the 
following bounds 


lla()I| < Boewv’ + & (14) 
eI < Bae + SE, (15) 
where the constants 6; > 0, i = 2,...,5 and vy, > v2 > 0 


are defined in the proof. 
Proof: It is straightforward to show that 7;(t) satisfies 
the dynamics 


Tui (t) + cmpi(t) + ypi(t)Hi(t) = —ypi(t)Zi+1 (8) 


—ypi(t)#;(t) + chy(t) + ha(t) (16) 


for all « = 1,...,n, where for the notational convenience 
we introduce a variable 7,4; = 0. The other variables are 
defined as p;(t) = y/ (t)y,(t), hi(t) = 8; (t)@, (t). Since 
a(t) and u(t) are bounded, it follows that «(t) is bounded. 
Therefore, there exist positive constants 6;, 62, 63 such 
that |lpi(Ilc. < 41, Wpilt)lle. < 52 and ||hi(t)\le., < 
63. Then, it can be concluded from equation (16) that c 
determines the damping and y determines the frequency 
of the signal 7;(t). It follows from the results of [5], that 
choosing c > 2,/5;7 damps the oscillations in 7(t) and 
guarantees the bound 
incall 
ial 
where 12 is proportional to ,/7 , and the positive constants 
B2, 64 and 65 are independent of y (see details in [5]). 
Substituting (9) we arrive to (14) with 63 = 6481 + 6365. 
Since A,. is Hurwitz, there exist positive constants dg and 
V2 such that the state transition matrix satisfies the inequality 
\|e4*|| < dge—%2". It follows that ||e(4-)4|| < dge— (tt, 
Then we obtain from (6) by the direct integration that 


Il*(4)|| < Bae~*" + d4]](4)|| + (17) 


\|z(t)|| < fs — lena = ge vate) 
tae [Miser ete Bae ake A , (18) 


since c and 1 are proportional to ,/7, which is much greater 
than v2 (adaptation is faster than the reference model). 


IV. CONTROL DESIGN 


In this section, we design four controllers, first three of 
which are used for the analysis purposes, and only the last 
one is implemented. 


A. Known System 


For the analysis purposes we first design a controller 
assuming that @ is known (ideal control). Following the 
standard backstepping procedure [4], we define new variables 


ze (t) =a ()-—of ,(t), i=1,...,n (19) 


(superscript 0 indicates that the variables corresponds to the 
ideal control signal) and stabilizing functions 


ag(t) =0, a(t) = —m(t) (20) 
a?(t) = —n(t) + a9_,(t), 1=2,...,n. 
The system (2) can be written in new variables as 
2°(t) = Az°(t) + blu®(t) — a2 (t)]. (21) 
Obviously, the control signal 
u(t) = kh! 2°(t) + ker (t) + an(t) (22) 


reduces the system into the reference model (3), hence the 
error e°(t) = 2°(t) — x,(t) satisfies dynamics 


e°(t) = A,e°(t). (23) 


Lemma 4.1: The controller defined by (19), (20) and (22) 
guarantees the control objective for the system (2). 

Proof: Since A,. is Hurwitz, it follows that e°(t) € Loo 
and e°(t) exponentially converges to zero. Hence, y°(t) = 
yr(t) + e9(t) exponentially converges to y,-(t). In addition, 
from r(t) € Lo it follows that x,(t) € Lo, which along 
with e°(t) € Lo. implies that z°(t) € £L... Boundedness of 
a(t), i=1,...,n and x°(t) is obtained recursively starting 
with x9 (t) = zP(t) € Loo. Then, u9(t) € Lo follows. 
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B. Certainty Equivalent Control 


Next, we design a controller for the identification model 
(4), by replacing the unknown parameter 6 with its estimate 
@(t) in the stabilizing functions 


Go(t) =0, dar(t) = —m(t) (24) 
dj (t) = —H,(t) + Qi_1(t), tS 250.15 
and introducing new variables as 


The identification model in new variables takes the form 


2(t) = AZ(t) + blu(t) — &n(t)] + c&(t). (26) 
Defining the control signal as 
a(t) = —k'2(t) + kpr(t) + Gn(t). (27) 
we obtain 
2(t) = A,2(t) + b,r(t) + cx(t), (28) 


which is in the form of modified reference model introduced 
in [5]. In this case the error signal é(t) = Z(t) — a,(t) 
evolves according to the dynamics 


é(t) = A,é(t) + c&(t). (29) 


Lemma 4.2: The controller defined by (24), (25) and (27) 
guarantees the control objective. 

Proof: According to Lemma 3.1, #(t) € Lo and 
z(t) > 0 as t > oo. Since A, is Hurwitz and r(t) € Loo, 
it follows that Z(t) € Lo, €(t) € Lo and é(t) > 0. 
Therefore, application of the controller defined by (24), (25) 
and (27) to both the system and the identification model 
results in 


y(t) = yr(t) + ert) + 1 (4) > y(t) 


In addition, from y;,.(t) € Loo it follows that 71(t) € Loo, 
hence @1(t) € Loo, since y, (a1) is continuous and 0; (t) € 
Loo according to Lemma 3.1. This implies that @2(t) € Lo 
and w2(t) € L5. Continuing this recursion we conclude 
that &(t) € Loo, @(t) € Lo, Gilt) € Lo, t= 1,...,n, 
ai(t) € Lo, t=1,...,n—1, and ti(t) € Lo. @ 

We notice that using the state x in the identification 
model (4) simplifies the stability analysis in the identification 
stage. However, the control design becomes problematic, 
because 4; (t) contains unknown parameter @ through the 
state derivative &(t). One way to overcome this issue is to re- 
place y;(a1,...,2;) with y, (#1, ...,%;) in the identification 
model (4), which brings additional terms j||@;||a ' (t)&(t) 
into the V(t) expression (8), where A; is the Lipschitz 
constant of y;(71,...,2;). In this case the same stability 
properties are guaranteed with the choice of c > 1+ 
ye Ail| 9: ||, and 4;(t) becomes implementable. However, 
the repetitive differentiations of 6;(t) introduces the multiple 
powers of the adaptation rate 7 into the control law, which 
makes the designer to keep the adaptation slow from the 
perspective of the control constraints. This is in conflict with 
the improvement of the performance by speeding up the 
adaptation. From this perspective, we use the command fil- 
tered backstepping approach from [2]. Although the method 
was introduce to simplify the process of determining the 
command derivatives in the backstepping procedure, it also 
allows to completely separate the identifier design from the 
controller design. Therefore, the identification process can 
be made as fast as the computational power allows. In the 
meantime, the high frequency oscillations associated with the 
fast adaptation are avoided with the proposed identification 
model by the proper choice of design parameters. 


(30) 


C. Command Filtering 


Following [2], we introduce the command filtered ap- 
proach for the design (19), (20) and (22), which will be 
used for the analysis purposes, and for the design (24), 
(25) and (27), which will actually be implemented. In the 
case of known system, the new variables are introduced as 
(superscript f indicates that the command filtered version) 


Hj Sa O02 440). tS 1yaains G1) 


where the command filter is designed as 
GP 1 (t) = wot a(t) 
Gf (t) = —2¢wo? »(t) a wlo? (4) = ad (t)] 
= Ope, 


(32) 


n—-1, 
with the initial conditions 0? , (0) = af (0) and o?,(0) = 0, 
and the stabilizing functions are defined as 
ag(t) =0, ay (t) = —m (4) 
ad (t) a —n} (t) se wat, 9(t), t= 2 
where the superscript f indicates that the corresponding 


quantities are computed when the command filtered control 
is in the loop. The system (2) is written in z-variable as 


(33) 


ee Oe 


2F(t) = Azf(t) + blu(t) — af (t)] + af (t) — @°(t), (G4) 
where 
al (t) 7,1 (t) 
al — : : o? — 
p at_,(t) a on—-11(t) 
0 0 


The compensated state is introduced as v°(t) = z/(t) — 
ip 
€°(t), where €°(t) is defined dynamically as 


O(t) =o )-a +80 (35) 


with: €°(0) = 0 for = 1,...,n— 1, and 620) = 0. In 
compensated state the system (2) takes the form 


#°(t) = Av°(t) + blu(t) — of (t)], (36) 
The control signal u(t) in this procedure is defined as 
uf (t) = —k' v(t) + k(t) taf (t). (37) 
Whereas the compensated error signal e2(t) = v°(t)—a,-(t) 
satisfies the dynamics 
ec(t) = Aree (t) (38) 


and obviously is exponentially stable, the uncompensated 
error ef (t) = zf(t) — w,(t) has dynamics 


ef (t) = A,ef(t) — af (t) + 6°(t). 


Lemma 4.3: The command filtered controller defined by 
(31), (32), (33), (35) and (37) guarantees the following 
relationships 


(39) 


ef (t)—e% (t)=O(e), €°(t) =O(e) (40) 
ofit)—of@=Ole), ¢=1,...,0-1 
wo? a(t) a(t) =O(e), +=1,. Lis 


where ¢ = 1/w (the proper choice of ¢ and w is discussed 
n [2]), and the notation O (¢) is adopted from [3] (p. 383). 
Proof: Although the error systems (38) and (39) are not 

in the standard backstepping form, the proof still follows the 
steps from [2]. | 
Now we design the last controller, which is actually im- 
plemented. It is the command filtered version of the certainty 


equivalent control, designed in the previous subsection. The 
uncompensated state is now introduced as 


zi(t) = 4 (t)-—oj-1i(t), t=1,...,n. (41) 
where o;_-1,1(¢) is the filter’s state given by 
di,1(t) = wo 2(t) (42) 
Gi,a(t) = —2woi,2(t) — wloia(t) — 4! ()] 
2=0,...,.n-1, 


with the initial conditions o;,;(0) = a! (0) and o;,2(0) = 0, 
and the stabilizing functions have the form 
ag(t)=0, at(t) =—Al 
éd (t) = —Af(t) toiie(t), i=2,...,n 
Here, we we introduce a short hand notation #/(t) = 
Tr 
6, (t)pf(t) and pf (t) = (xf (t),..., x7 (t)). The iden- 
tification model in z-variables takes the form 
z(t) = Az(t) + blu(t) — ah (0) 
+cH(t) + f(t) — a(t). (44) 
The compensated state and its dynamics for the identifi- 
cation model are similarly introduced 


E(t) = o1a(t) — Af (t) + G41 (4) (45) 
with €;(0) = 0 fori =1,...,n—1, &,(t) = 0, and v(t) = 


(43) 


2(t) — €(t), 
0(t) = Av(t) + blu(t) — af (t)] + cx(t). (46) 
The control signal to be implemented has the form 
u(t) = —k' v(t) + kpr(t) + df(t). (47) 
The resulting compensated error signal e,(t) = v(t) — x,-(t) 


and uncompensated error signal e(t) 
dynamics 


= z(t) — x,(t) satisfy 


€.(t) = A,e.(t) + cx(t) (48) 
é(t) = A,e(t) — &(t) + &(t) + cx(t). (49) 
Lemma 4.4: The command filtered controller defined by 


(41), (42), (43), (45) and (47) guarantees the following 
relationships 


e(t)—e(t)=Ole), &(t) =O(e) (50) 
oii(t) — di(t) =O(e), t=1,...,n-1 
woio(t) — di(t)=O(e), i=1,....n-1. 


Proof: Since the exponential convergence of &(t) is 
not guaranteed, Tikhonov’s theorem ( [3], Theorem 11.2) 
cannot be directly applied to the system comprised of (44), 
(45) and (49). However, since «(t) does not depend on «, a 
simple state transformation s(t) = e(t)— p(t), where ju(t) is 
dynamically defined as f(t) = A,ps(t) +cx(t), results in the 
system 8(t) = A,.s(t)—@(t)+@(t), for which the hypothesis 
of the Thikhonov’s theorem can be verified following the 
steps from [2], and the last two relationships in (50) can be 
concluded along with s(t) — s(t) = O(c), where s(t) = 
é(t) — p(t) and satisfies the exponentially stable dynamics 
s(t) = A,8(t). It follows that e(t) — é(t) = s(t) — 8(t) = 
O (e), which completes the proof. a 


V. PERFORMANCE ANALYSIS 


The following two lemmas are needed to prove our main 
result. 

Lemma 5.1: Let the command filtered controller for sys- 
tem (2) be defined by (31), (32), (33), (35) and (37). Then 
all closed-loop signals are bounded and 


at (t) — #°(t)=O(e). (51) 
In addition, if w is sufficiently large, then 
uf (t) — u(t) = Ole) . (52) 


Proof: Since e°(t) € Loo it follows from (40) that 
ef(t) € Loo, implying that z/(t) € Lo. It follows form 
Lemma 4.1 and Lemma 4.3 that a? (t) € Lo and of 4(t) € 
Lo for i = 1,...,n —1. Then, (31) implies that a(t) € 
Loo. Therefore nf(t) € Loo, hence ad (t) € Loo for all 
i =1,...,n. Since e,(t) is exponentially stable, it follows 
that v°(t) € Lo, and hence u(t) € Lo 

It is straightforward to compute the difference 


wf (t) — w(t) = ef (t) — eP(t) + 09s (t) — a(t). 


Since ef (t) — ef(t) = O(e) and of (t) — aP(é) = Oe), 
it follows that a! (t) — x9(t) = O(c) for all i =1,...,n 
Next, we compute the difference 

uf (t) — u(t) = —k" [w(t) — 2°(8)] + af (8) — a(t). 


Since €°(t) = O (e), it follows that v°(t) — z°(t) = O (e). 
On the other hand af(t) — a®(t) = nf(t) — m(t) + 


woy_19(t) — a? _4(t). Since ~,,(@) is smooth, we have 
vy, (at (t)) — y,(x°(t)) = O(c). Then it follows that 
ut (t) — u°(t) = O (e) if w > max(||kl|, |]On/I). | 


Lemma 5.2: Let the command filtered controller for sys- 
tem (2) and identification model (4) be defined by (31), (32), 
(33), (35) and (37). Then all signals are bounded and 


&!(t)— #(t) =O(e). (53) 
In addition, if w is sufficiently large, then 
ait (t) — G(t) = O(e). (54) 


The lemma is proved similar to the previous one. 

Theorem 5.1: Let the system’s controller be defined ac- 
cording to command filtered scheme given by (41), (42), 
(43), (45) and (47). Then the input and output tracking errors 
satisfy the following upper bounds 


|d(¢)| < Boe" + + Oe) (55) 
lel) < Bee"* + B+ Ole), (56) 


where ti(t) = ||u°(t) — u(t)||, 86, G7, Bs, Go and v are 
positive constants defined in the proof. 
Proof: It is easy to see that 


u(t) = u(t) — uf (t) + ul (t) — a(t) 

= O(c) — kT a(t) + af (2), 
where w(t) = v°(t) — v(t) and a, (t) = af(t) — af(t). 
Obviously, (t) satisfies the dynamics 
v(t) = A,o(t) + cx(t) 


with the initial conditions v(0) = (0) — (0), where we 
denote 7;(t) = ni (t) — fy (t) = 0;(t)p! (t) for each i = 
1,...,n. similar to (18), one can obtain form (58) that 


sae —vVat Bir 
I|o(¢)| < Bioe ze Vy? (59) 
Systems response 
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Fig. 1. Adaptive step response and corresponding adaptive control signal 
vs ideal step response and ideal control signal. 


Next, we observe that the signals 4;,1(t) = o?  (t)—o;,1(t) 
and G;,2(t) = o)5(t) — 0i,2(t) satisfy the operator equations 


Gia(s) =Gi(s)él(s), Gils)=apeacr (60) 
Gi,2(s) = Go(s)af (s), Ga(s) = FT atestuF ? 
where @/ (t) = af (t) — af (t). Since ||Gi(s)||42,, = 1 for 


! it follows from [3] 


¢ > V2/2 and ||G2(s) lla. 
(p. 201) that 


= (2¢)~ 


Fallen < laf Olle. 
F.2lle. < 2A Ol cu 
Now we can recursively compute the bound on a/(t) 


using the definitions (33) and (43). For 7 = 1, we have 
a! (t) = al (t), therefore 


(61) 


Bs 


|@i(t)| < Boe" + Vi (62) 
For i = 2,...,n, we have &! (t) = —f;(t) + w6i,2(t), hence 
in (| < anBaem + (63) 


n-1 
where gq, = 1+ aor ~+ i) 
Combining the relationships (57), (59) and (63) we obtain 


(55), where 86 = ||K||F10+¢nB2, 87 = ||| G11 +9n 8s, since 
V, > V2 for large values of yy (fast adaptation). 
To prove (56), we notice that y(t) — y,(t) = oy a 


a! (t)— 2,1 (t) = #f (t) + e1(t). Since ey (t) = a(t) + O(e), 
it follows that y(t) — y,(t) = Z(t) + é1(t) + O (e). Using 
(18) one can obtain from (29) that 

EWI < Brze* + 8 (64) 


Vy? 


and the relationship (56) follows. The proof is complete. 

It follows form Theorem 5.1 that the bounds on the input 
and output tracking errors can be systematically decreased 
by choosing large values for w and +. 
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Fig. 2. Estimates of 7 and 72 vs true values for the step input. 
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Fig. 3. Adaptive sinusoidal response and corresponding adaptive control 
signal vs ideal step response and ideal control signal. 


Remark 5.1: The proposed approach can be readily ex- 
tended to multi-input-multi-output systems with unmatched 
uncertainties. Also, time variant parameters and external 
disturbances can be introduced in the proposed approach. 


VI. SIMULATION RESULTS 


As a simulation example we consider a third order system 


&1(t) = ro(t) + 0127 (t) (65) 
&2(t) = x3(t) + O2x1 (t)x2(t)) 
&3(t) = u(t) 
with the unknown parameters 6; = 3 and 62 = 2. The 


reference model is selected with k' = [8 10.4 5.2] and 


k,, = 8. The identification model is designed with c = 495 
and yy = 500000. For the command filtering we set w = 500 
and ¢ = 0.8. Figure | displays systems response to step input 
along with the command filtered certainty equivalent control 
signal vs the systems step response (ideal response) and 
control signal (ideal control) for the conventional backstep- 
ping design assuming parameters are known. The identifier’s 
performance is displayed in Figure 2. It can be observed 
that good tracking is achieved in all signals for the unit 
step command. The proposed controller achieves a similar 
performance for the sinusoidal command as can be seen from 
Figures 3 and 4. 
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Fig. 4. Estimates of 71 and 12 vs true values for the sinusoidal input. 


VII. CONCLUDING REMARKS 


We have presented an indirect adaptive control method for 
nonlinear systems with unmatched uncertainties that follows 
the certainty equivalence principle. The approach uses a fast 
identification model, which is independent of the control de- 
sign and achieves desired transient and steady state properties 
by the proper choice of the design parameters. The controller 
is in the form of the command filtered backstepping control. 
The resulting tracking errors can be decreased as desired by 
speeding up the adaptation and command filtering processes, 
subject to available computational power. 
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